Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 13 September 2011 (MN WTp^ style file v2.2) 



Modelling the VHE flare of 3C 279 using one zone leptonic 
model 



o 

(N 



S. Sahayanathan^* and S. Godambe^f 

^ Astrophysical Sciences Division, Bhabha Atomic Research Centre, Mumbai - 400085, India 
^Department of Physics and Astronomy, The University of Utah, Salt Lake City, Utah-84112, 



USA. 



13 September 2011 



X 
6 



> 

(N 
(N 

o 



ABSTRACT 

We model the simultaneous observations of the flat spectrum radio quasar 3C 279 at 
radio, optical, X-ray and very high energy (VHE) gamma ray energies during 2006 
flare using a simple one zone leptonic model. We consider synchrotron emission due to 
cooling of a non-thermal electron distribution in an equipartition magnetic field and 
inverse Compton emission due to the scattering off synchrotron photons (SSC) and 
external soft photons (EC) by the same distribution of electrons. We show that the 
VHE gamma ray flux cannot be explained by SSC process thereby suggesting the EC 
mechanism as a plausible emission process at this energy. The EC scattering of BLR 
photons to VHE energies will be in Klein-Nishina regime predicting a steep spectrum 
which is contrary to the observations. However the infrared photons from the dusty 
torus can be boosted to VHE energies with the scattering process remaining in the 
Thomson regime. Though EC process can successfully explain the observed VHE flux, 
it require a magnetic field much lower than the equipartition value to reproduce the 
observed X-ray flux. Hence we attribute the X-ray emission due to SSC process. We 
derive the physical parameters of 3C 279 considering the above mentioned emission 
processes. In addition we assume the size of emission region constrained by a variability 
timescale of one day. This model can successfully reproduce the broadband spectrum 
of 3C 279 but predicts substantially large flux at MeV-GeV energies. 
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1 INTRODUCTION 

Flat Spectrum Radio Quasars (FSRQ) are radio loud 
Active Galactic Nuclei (AGN) with broad emission lines 
and a non-thermal spectrum extending from radio to 
gamma ray energies. They are characterized by luminous 
core, rapidly variable non-thermal emission, high radio and 
optical polarization, flat-spectrum radio emission and/or 
superluminal motion. Similar properties are also observed 
in BL Lacs and these two types of AGN are classified as 
blazars. According to AGN unification scheme, blazars are 
the class of AGN with a rel ativistic jet pointed close to the 
line of sight of the observer (|Urrv fc Padovanil (lOoF) ) . The 
spectral energy distribution (SED) of FSRQ is characterized 
by two broad humps. The low energy hump in the SED is 
well understood as synchrotron emission from a relativistic 
distribution of electrons. Whereas, the origin of high energy 
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hump is still matter of debate. There are several models 
to explain the high energy emi ssion based on either lep- 
tonic o r hadronic interact i ons IIBlandford fc Levinson. A.I 
("logs'): 'Bloom fc Marscher' ("iQge'): 'Aharonian et al.1 (|200of )7 
Pohl, & Schhckciscr (2000); Mannheim (1998)). Under the 
assumption of leptonic models, high energy emission is 
explained via inverse Compton scattering of soft target 
photons. The target photons can be t he synchrotro n 
photons (synchrotron self Comp ton (SSC)) ('Konigl' (ToSl]); 
Mar scher fc Gear (.1985. ): ,Ghisellini fc Mar aschi (198E1)) 
the ph otons external to th e jet ( external Comp- 
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ton (EC ))(lBegelman fc Sikoral (Il987l): . 
(jl989l ): iDermer. Schlickeiser. fc Mastichiadi? 
These external phot ons can be the accre - 
tion disk photons (Dermer fc Schlickeiserl (|l993l ): 
iBoettcher. Mause. fc Schlickeiser! ( 19971 )) or the accretion 
disk photons reprocessed by b road fi ne region (BLR) clouds 
(Sik ora. Begelman. fc Reej (Il994l ): iGhiseUini fc Madaul 
( 199( |)) or the infrared r adia tion from the dusty torus 
(jSikora, Begelman. fc Ree3 (| 19941 ); iBlaieiowski et all 
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Ghisellini fc Tavecchid (|2008l ')'). Hadronic models 
explain the high energy emission as an outcome of the 
synchrotron proton emission and proton-photon interac- 
tions with the synchrotron photo ns (synchrotron proto n 
blazar model) jMannheiml ([l998l); iMticke et all l|2003h : 
iBottcher. Reimer. fc Marscheil (|2009l )). 

3C 279 {z = 0.536) is a well studied source at var- 
ious energy bands with several s imultaneous multi wave- 
length campaigns(|Maraschi et aP (fl994): ,Haxtman et"al] 
(| 19961 ): IWehrle et all (|l998l ): iBdttcher et all (|200'it l). It 
was the first blazar observed at gamma ray energies by 
the satellite based experiment EGRETfEnerg e tic G amma- 
Ray Experiment Telescope) (|Hartman et al.l (| 19921 )) and 
the first FSRQ to be detected at GeV-TeV (very high 
energy (VHE)) gamma ray energies by ground based 
atmospheric Cherenkov experiment MAGIC (M a-jor At- 
mosph eric Gamma-Ray Imaging Cherenkov ){ Alb ert et alj 
l|2008h ). The flux of 3C 279 is also known to be strongly 
variable at ra dio, infrared, o ptical , UV, X-ray and gamma- 
ray energies l|Makino et alJ l|l989l l and references therein). 
During 2006, 3C 279 undergone a dramatic flare and 
was observ ed by WEBT ( Whole Earth blazar Telescope) 
campaign jBottcher et all (|2007l )) at radio, near infrared 
and optical frequencies and simultaneously monitored by 
Rossi X-ray Tim i ng E xplorer (RXTE) at X-ray energies 
IChatteriee et all (120081)) a nd by MA G/C telescope at VHE 
gamma ravs (jAlbert et al.l (j2008i )). 

The SED of 3C 279 is modelled as synchrotron and 
inverse Compton emission using these simultaneous obser- 
vations at various energies. The explanation of the Compton 
hump as an outcome of SSC process due to the cooling of 
a relativisti c power-law electron distribut io n is not very 
successful (jBottcher. Reimer. fc Marsche^ ()2009l )). This 
interpretation had a serious drawbacks in explaining the 
high energy gamma ray (> IGeV) detection from 3C279 
since it require very high rad iation energy density compared 
to ma gnetic energy densitv (|Maraschi. Ghisellini. fc Celottil 
j 19921 )). It also require unus ually low magnetic fleld to ex- 
plain the VHE emission l|Bottcher. Reimer. fc Marscheil 
l|2009l )). However in case of BL Lac objects, 
the models considering the SSC process are able 
to reproduce gamma ray e m ission successfully 
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1 19991 ): iBhattacharvva. Sahayanathan. fc BhattI l|2005l )). 
Other alternative explanation for the Compton hump 
of 3C 279 is the EC mechanism with the target photons 
external to the jet. Target photons from the accretion disk 
will be strongly deboosted due to high Lorentz factor of 
the jet (r ~ 10 — 25) and may not play an important 
role (however see iBoettcher. Mause. fc Schlickeisen (| 19971 ) 
for an alternative). The photons from BLR clouds cannot 
be boosted to VHE gamma rays since the scattering 
of Lyman alpha li ne emission (t he d ominant emission 
from BLR clouds (iFrancis et~al] (|l99ll))) will result in 
Klein-Nishina regime (|Ghisellini fc Tavecchid (|2009l ): sec- 
tion i|3.2p predicting a steep spectrum contrary to the 
observed hard spectrum. The third option, infrared photons 
from the dusty torus, can be a plausible candidate for 
the target photons in EC process since the scattering of 
these photons to V HE gamma ray energie s will still be in 
Thomson regime l|Ghisellini fc Tavecchid (|2009l ): section 
i|3.2|) . The importance of Compton scattering of infrared 



photons from du st y tor us was first studied in detail by 
iBlazeiowski et all l|2000l ) to explain high energy gamma 
rays detected by EGRET. Though EC process is widely 
accepted to explain the g amma ray emission of 3C 279, 
iLindfors. Valtaoia. Tiirle^ (|2005l ) argues in favour of SSC 
mechanism based on observed radio-gamma ray correlation 
and the quadratic dependence of the synchrotron and 
inverse Compton peak fiuxes during the flare. 

In the present paper we use the simultaneous obser- 
vation of 3C 279 at radio, near infrared, opt ical. X-ray and 
very high ener gy (VHE) gamma ray energies JBottcher et al.l 
(2007); Chatte riee et al.l (|2008l ): [Albert et al.l (|2008l )) to de- 
duce the physical parameters of the source. We consider 
synchrotron, synchrotron self Compton (SSC) and external 
Compton (EC) emission from a power-law distribution of 
electrons. We assume a magnetic field which is in equiparti- 
tion with the particle energy. In the next section we outline 
the model and in section iJ3]we show the plausible target pho- 
tons for the inverse Compton process to explain the observed 
flux at different energy bands. In section ^we present the 
estimated physical parameters of 3C 279 and discuss impli- 
cations of the present model and the results. A cosmology 
with = 0.3, f^A = 0.7 and Ho = 70 km Mpc'^ is 
used in this work. 



2 MODEL 

We assume the emission region to be a spherical blob moving 
down the jet at relativistic speed with Lorentz factor F. Since 
the jets of blazars are aligned towards the line of sight of 
the observer, we assume the Doppler factor 5 ~ F. The 
radiation from emission region is due to synchrotron and 
inverse Compton cooling of relativistic electrons described 
by a broken power-law distribution (quantities with prime 
are measured at the rest frame of the emission region) 
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where 7' is the Lorentz factor of the electron, p and q are the 
power-law indices before and after the break corresponding 
to the Lorentz factor 7^ and K is the particle normalisation. 
In the blob frame protons are assumed to be cold and con- 
tribute only to the inertia of the jet. The emission region 
is permeated with tangled magnetic fleld B. If we assume 
the equipartition condition between the magnetic fleld en- 
ergy density((7B) and the particle energy density, then the 
equipartition magnetic field Beq will satisfy 
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where m is the electron rest mass and c is the velocity of 
light. We neglect the contribution of protons in equation ((2]) 
due to the assumption they are cold and do not contribute 
to the particle energy density. The size of the emission re- 
gion R' can be approximated from the observed variability 
timescale tvar as 

R ~ ; rCtvar (3) 

(1 + ^) 

where z is the redshift of the source. 
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3 INVERSE COMPTON PROCESS IN 3C 279 

The radio-to-optical emission of 3C 279 is considered to be 
synchrotron emission and X-ray-to-7-ray emission is pro- 
duced by inverse Compton scattering of soft photons. The 
soft target photons can be the synchrotron photon itself 
(SSC) and/or the photon field external to the jet (EC). 



3.1 Synchrotron Self Compton (SSC) 

In SSC mechanism, the synchrotron photons are boosted to 
high energies via inverse Compton scattering by the same 
electron distribution responsible for the synchrotron emis- 
sion itself. An approximate expression for synchrotron emis- 
sivity e'sy„{i'')[ergs cm~^ Sr~^ H z~^] due to an electron 
distribution N'{'y') can be obtained by describing the single 
parti c le em issivity in terms of a delta function 6(1^' — 7'^«^l) 
Isiiul (|l99ll ')') (Appendix El) 
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Where ctt is Thomson cross section, vl = tt^— is 

' 27rmc 

the Larmor frequency and v' is the frequency of the 
emitted photon in the rest frame of the emission re- 
gio n. The observed flux Fayn(v) \ erqs cm~'^ H z^^\ will 
be (|Begelman. Blandford. fc Reed l|l984 )') 
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Where V' is the volume of the emission region and dL is the 
luminosity distance. The rising part of the observed syn- 
chrotron flux in the SED due to the particle distribution 
given by equationU]) can then be written as 

Fsyniu) ~ s(z,p)S'^ R'^Ku-"^ Jy (6) 

for ly < Sj'i^VL/ii + z). Here s{z,p) is a function of p and 
z. For z = 0.536 and p — 2.02 (corresponding to a radio 
spectral index 0.51), s = 1.7 x 10"^^. 

The observed synchrotron and SSC frequency corre- 
sponding to the peak flux in the SED can be approximated 
as 
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Using equations©, Q, Q and (|9j (for p = 2.02) we get 
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€'aa^{iy')[ergs cm~^s~^ Sr~^ Hz~^] for the particle dis- 
tribution given by equationQ can be approximated as 
(Appendix |B| 
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where O is the Heaviside function and 
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However the gamma-ray peak frequency expressed in equa- 
tion (|10p is too low to interpret the observed VHE flux as an 
outcome of SSC mechanism (see Figll}). 

The emissivity due to SSC mechanism 



The observed SSC flux can be obtained following equa- 
tion(lS| by replacing e'gy„{u') with e!,sc('^')- H should be 
noted here in deriving equations H12|) - p6[) we have as- 
sumed the scattering in Thomson regime. However at 
VHE energies this may not be valid since Klein-Nishina 
effect will be more prominent and sh ould be included 
(|Tavecchio. Maraschi. fc Ghisellinil (|l998l ')'). 



3.2 External Compton (EC) 

In EC mechanism the target photons for the inverse Comp- 
ton scattering are accretion disk photons, reprocessed disk 
photons by the BLR clouds and infrared photons from the 
dusty torus. In the following discussion we consider photons 
from the BLR clouds and the dusty torus only, since the 
photons from the accretion disk will be strongly redshifted. 

If the VHE emission is produced due to inverse Comp- 
ton scattering happening in Klein-Nishina regime, then the 
VHE spectrum should be steeper than the observed syn- 
chrotron spectrum. However the intrinsic VHE spectrum 
of 3C 279 is substantially harder after correction for the 
absorpt ion effects due to extra galactic background light 
(EBL) (lAlbert et all (120081)'). CThough the co nsidered EBL 
model JPrimack . Bulloc k fc Somervilld (|2005l ')') is debatable 
(|Stecker fc Scully, (,2009i' )'). Therefore, from the inferred in- 
trinsic VHE spectrum, one can conclude that the scattering 
process must be in Thomson regime. 
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iGhisellini fc Tavecchid (|2009l ) showed that the scatter- 
ing of BLR photons to VHE energies will be in Klein-Nishina 
regime. The condition for Thomson scattering is given by 

^ < 1 (17) 

Where h is the Planck constant and is the frequency 
of target photon in AGN frame (quantities with star are 
measured at the AGN frame) . The frequency of the scattered 
photon in the observer's frame is 



From equations (|17|) and (|18|) we get 
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The dominant emission from BLR region is th e line emission 
at 2.4 7 X 10^^ Hz corresponding to Lyman-a (|Francis et al.l 
l|l99ll )') and hence the scattering will be in Klein-Nishina 
regime. However the scattering of infrared photon (10^^ — 
lO^* Hz) from the dusty torus to VHE energies can still re- 
main in Thomson regime. The interpretation of VHE emis- 
sion as a result of inverse Compton scattering of infrared 
photons is also used to explain the VHE emission from 
the intermediate BL La c objects, W C omae and 3C 66A, 
detect ed by VERITAS (|Acciari et al.l (|2008l ): lAbdo et all 

dmS)). 

The EC emissivity €'^c{u')[ergs cm ^Sr ^Hz ^] 
due to scattering of an isotropic mono energetic photon field 
of frequency v* a nd energy dens ity u* (in lab frame) can be 



or irequency u a nd energy aens ity 
approximated as jPerme jf (|l995l ')') 
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The observed EC flux can be obtained following equation([5]) 
by replacing e'^y^{v') with e'^^iv'). 

If the observed Compton flux can be represented by a 
broken power-law 
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with a = ^-Y- and /3 = then the peak Compton fre- 
quency in SED will be 



(22) 



where vi < fp,ec and V2 > i'p,ec- Since Vp^^c corresponds to 
the scattering off the target photon by the particle of energy 
'ylmc? , we can write 

,2 * 

t'p.ec ~ ^ ^-^ Ib V (23) 

If we assume both the X-ray and the VHE emission are 
produced by the EC scattering off infrared photons itself, 
then using equations ((Tj, (|22|l and (|23p . we can write 

,26, 



B ~ 5.4 X 10 



5 X 10" j \2.2 X 10" 
• (Fee(i^i)/1.5 X 10-''jy)(t.i/10^^) 
XF^c(y2)l\n X 10-12 Jy)(!/2/2.1 X 1025) 



G (24) 



Here we have used q = 0.51 and /3 = 1.6 corresponding to 
observed spectral indices at radio/X-ray and optical ener- 
gies. However the equipartition magnetic fleld deduced from 
equations((2]), ([3]) and ((6| is 
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Hence we require a magnetic fleld much lower than its 
equipartition value to reproduce the observed X-ray flux 
through EC scattering of the infrared photons. However if 
X-ray emission is due to SSC process, the required magnetic 
field will be close to its equipartition value. 

Based on these arguments, we model the broadband 
spectrum of 3C 279 considering synchrotron, SSC and EC 
processes. The synchrotron emission is dominant at ra- 
dio/optical energies whereas SSC is dominant at X-ray en- 
ergy and EC at VHE gamma ray energy. 



4 RESULT AND DISCUSSION 

The main parameters governing the observed broadband 
spectrum of 3C 279 are the B, K, p, q, 7^, R' , T, v* and it*. 
Out of these the indices p and q can be constrained using the 
observed spectral indices at radio and optical energies. If we 
assume the target photon energy density u* from the dusty 
torus as a blackbody, then the peak target photon frequency 
V* can be written in terms of u* as 



= 2.82 



Kb 



(26) 



where Kb is the Boltzmann constant and ct is the Stefan- 
Boltzmann constant. The remaining six parameters B, K, 
75, R! , r and u* can be calculated using equations®, @, 
(0,0, (fTTj) a nd (IM . We consider o ne day as the variability 
timescale t^ar (jHartman et al.l (|l996l )). The derived parame- 
ters are given in table Q and the corresponding broadband 
spectrum is shown in figure ([T|. We have chosen 7m!n = 40 
and 7mai: = 1.0 x 10^ to extend the observed spectrum from 
radio to VHE energies. The approximate expression for syn- 
chrotron, SSC and EC fluxes (discussed in earlier section and 
AppendixS&lBj are used to derive the parameters whereas 
for reproducing the observed flux in figure ([1]) the exact ex- 
pressions are used and solved numerically. With the given 
set of parameters if we consider the number of cold protons 
are equal to the number of non-thermal electrons, then the 
energy density due to protons will be ~ 1.5 x 10~^ ergs / cm? . 
This is an order less than the electron energy density (table 
0) and hence its exclusion will not affect Beq considerably. 

From the figure one can find that model predicts 
excessive fiux at MeV energies which is nearly an order 
more than the highest flux detected firom 3C 279 by EGRET 
during its entire mission (|Hartman et all (|200lh '). The flux 
obtained during 2008 observations by Fermi is also much 
lower. Hence it can be argued that this prediction may be 
unlikely. However there were no simultaneous observations 
of source at MeV-GeV energies during the WEBT campaign 
and MAGIC detection and hence such a dramatic flare can 
not be ruled out. 

The deduced temperature of the dusty torus (T = 
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Table 1. Model Parameters 



Parameter 


Symbol 


Numerical Value 


Particle spectral index 


P 


2.02 


(low energy) 






Particle spectral index 


g 


4.2 


(high energy) 






Magnetic Field 


Beq 


0.67 G 


( equipartition ) 






Particle energy density 




1.8 X 10~^ ergs/cm^ 


Particle spectrum break energy 




832 


(units of mc^) 




Emission region size 




4 X lO^S cm 


Bulk Lorcntz factor 


r 


26 


IR dust temperature 


T* 


860 K 


IR dust photon energy density 


u* 


4 X 10~^ ergs/cm^ 




Log(v)(Hz) 

Figure 1. Spectral fit to the SED of 3C 279. The radio(plus), op- 
tical{cross) and iJXTE X-ray data (butterfiy) are obtained from 
Boettcher et al. (2007) and the VHE data is reproduced from Al- 
bert et al. (2008). The Chandra (filled circle) and Swift (open cir- 
cle) X-ray data and the INTEGRAL (triangle) data are obtained 
from CoUmar et al. (2010). The MeV-GeV data represented by 
open boxes are the iJGRET data obtained from Hartman et al. 
(2001). The MeV-GeV data represented by filled boxes are _Fermi 
data obtained during 2008 (Abdo et al. (2010)). The dashed curve 
is the synchrotron spectrum, the dotted curve is the SSC spec- 
trum and the dot-dashed curve is the EC spectrum. The solid 
curve is the total contribution of the different emission processes. 



860A") is consistent with the temperature obtained from 
the inner region of the dusty t orus (> 800K) of the nearby 
Seyfert type 2 AGN NGC 1068 (|jafre et al.l l|2004l )1. It is also 
close to the temperat ure range (900 — 1300K) suggested by 
Ide Vries et al] (|l998l ) based on near-infrared modelling of 
gigahertz peaked spectrum (GPS), compact steep spectrum 
(CSS) and Fanaroff & Riley type II (FR II) sources. 

The emitted VHE photons can be absorbed by soft tar- 
get photons through photon-photon pair production mech- 
anism. The condition on photon energies to be transparent 



against pair formation is 

{1 + z){hu^){hiy*) < {mc^f 



(27) 



If vvHE ~ 10 Hz then it can pair produce with tar- 
get photons with frequency ly* > lO^^i^z. Hence the 
emitted VHE photons can pair produce with the BLR 
photons, however the infra red photons from the dusty 
torus will be below the threshold. If the emission re- 
gion is within the BLR region, the VHE photons can 
be a bsorbed by the BLR photons throu g h pair produc- 
tion (Bottcher. Reimer. fc Marscheil l|2009l ): iBai. Liu.fc Mai 
l 200m . Hence the detection of 3C279 at VHE gamma 
ray energies demand the emission region may be ahead 
of the BLR r e gion to avoid the severe 77 absorption 
iBai. Liu.fc Mai 



severe 

absorption of VHE 



(|2009l ) studied the 77 
gamma rays from 3C 279 and suggested the emission region 
must be located within the BLR region. However the present 
work (and also Bottcher, Reim er. fc MarscheJ (2009)) re- 
quire the emission region to be far from the BLR region to 
reproduce the simultaneous broadband spectrum of 3C 279. 

If we consider the dusty torus as a n ann ular ring cover- 
ing the central source (|Pier fc KrolikI (|l992l )) then the dis- 
tance of the inner wall of the torus (Rir) can be estimated 
as 



RiR ~ 



1.6 



(- 



1/2 



parsec (28) 



where L„„ is the accretion disk luminosity. Also, if the dis- 
tance between BLR clouds and the central so urce (Rblr) 
of 3C 279 satisfy the size- luminosity relation of lKaspi et al.l 
(l2007l l. then 



Rblr ~ 0.052 



/ ALa(1350 A) 
V lO^B erps s-i 



parsec 



(29) 



For the considered variability timescale of 1 day 
(jHartman et al.l (|l996l )'). the distance of the emission region 
from the central source (Rfi) can be approximated as 



Rfl ~ CtuarF 



parsec 



(30) 



which is ~ 10 times farther than Rblr supporting our infer- 
ence. Also from eauation (|28p the emission region lies within 

the dusty torus i.e Rblr < Rfi < R iR- 

iBottcher. Reimer. fc Marscheil ([2009') suggested a 
multi-zone leptonic model as a possibility to explain the 
VHE emission from 3C 279 since one zone model require 
low magnetic fields and/or very high Lorentz factor of the 
jet. However they arrived to this conclusion by considering 
the BLR photons as the target photons for the external 
Comp ton process instead of inf rared photons from the dusty 
torus. iBla^eiowski et all (|2000l ) considered a model which is 
similar to the one described in the present work, but their 
aim was to project the importance of the Comptonisation 
of infrared photons from the torus to reproduce the high 
energy spectrum. In the present work, we have used 
simultaneous observation of 3C 279 at different energies to 
deduce the physical parameters of the source. 
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5 CONCLUSION 

We reproduce the observed simultaneous broadband spec- 
trum of 3C 279 using a simple one zone leptonic model con- 
sidering synchrotron, SSC and EC processes. From the ra- 
dio/optical synchrotron spectrum we show that the VHE 
emission cannot be attributed to SSC process whereas EC 
scattering of IR photons from the dusty torus can explain 
the observed VHE emission. Interpreting the X-ray emission 
as continuation of EC spectrum require magnetic field much 
lower than its equipartition value. However an explanation 
based on SSC origin of X-ray require a magnetic field which 
is comparable to the equipartition value. 

The model predicts large fiux at MeV-GeV energies. 
Since the data in these energy range is not available during 
the considered flare, such a prediction cannot be ruled out. 
The model parameters describing the source are estimated 
by considering a magnetic field which is in equipartition with 
the particle energy density. A deviation from this condition 
will reflect in the value of the estimated parameters. Also 
if the jet matter contains considerable amount of energetic 
protons, then the contribution from these protons should 
be included in equipartition magnetic field. However in the 
present work, the proton contribution to the total particle 
energy is considered to be negligible. Also the size of the 
emission region estimated from equation(12J is only an upper 
limit and the actual size may be smaller. Moreover in reality 
the size of the emission region may be different for different 
energy bands. These variation in the emission region size 
will also be reflected in the estimated parameters. 
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APPENDIX A: SYNCHROTRON EMISSIVITY 



The SSC emissivity will then be 



The synchrotron emissivity esyn{v) at frequency u due to a 
relativistic electron distribution A^(7) can be written as 



(Al) 



Where Psyn{j,i^) is the average synchrotron power emitted 
by an elec tron with energy 'ymc? at frequency v. Following 
IShul (|l99ll ). we can express Psynij, v) as 



Psyn{n,v) = g/3^7^CCrT[/s</>^(7) 



(A2) 



where /3c(« c) is the velocity of the electron and the function 
4>v{')) satisfies 



/ 



If we approximate (j>v('y) as a delta function 



(A3) 



(A4) 



Then we can perform the integration (IA1[) and the syn- 
chrotron emissivity will be 



487r 



(A5) 



APPENDIX B: SSC EMISSIVITY 

The SSC emissivity essc(/^) at frequency u due to a relativis- 
tic electron distribution A''(7) can be written as 



£ssc ( 



oo 



Pssc{l,p)N('y)d'y 



(Bl) 



Where Pssc(7, i^) is the average SSC power emitted by an 
electron with energy ^m(? at frequency v and can be written 
as 



Pss.(n,v) = -^7 car / V4^,7) (B2) 



^ssA^) ~ ^coT ju(^^] iV(7)d7 (B6) 
1 



If we write 



.(0 



(B7) 



then from equation (|A5|l we get 

Cx 

Pc 2 2 ~ — — 

^ ' 367r2 ^ ^ 



-^f-\/?)^W (B8) 
7 \1\ vl) 



For a broken power law electron distribution (equation([T| 
we get 

Rc 



where 



^^^^K'cylB'v^-^v-^jiv) (B9) 



V72/ 9-P V^^i 



76 V 



8(72 — 71) 



+ 



In 



74 \ , 76 / 



73 / q-p \VL 



^174 Imax I 7™aa; 1/ 



here is the Heaviside function and 



71 = MAX 7„i„, — 



72 = M7iV 7b 



73 = MAX 76 



72 = M7Af -fma: 



7i> V '^i- 



8(74 - 73) 

(BIO) 

(Bll) 
(B12) 
(B13) 
(B14) 



where 



00 



(B3) 



is the energy density of the target photons (in this case syn- 
chrotron photons) and ^ is the frequency of the target pho- 
ton. The function tpi^i^j'y) will satisfy the condition 



7/'^(C,7)di^ = 1 



(B4) 



Since the frequency(i^) of the photon scattered at Thomson 
regime is ~ 7^^, we can approximate ^i/(C)7) as a delta 
function 



(B5) 



